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We present results on the spectrum of B and D mesons including radial excitations and discuss the pseudoscalar 
decay constant. The results are obtained at /3 = 5.7 in the quenched approximation using NRQCD for the heavy 
quark. To study scaling violations we also compare to results obtained at j3 = 6.0. 



1. INTRODUCTION 

We want to detail our new and as yet pre- 
liminary simulation results on the spectrum and 
the renormalised decay constant of heavy light 
mesons. They were obtained on coarser lattices 
compared to previous analyses, see Q for a re- 
view. This allows us to study scaling violations. 
A second motivation is the possibility of reaching 
the charm region with NRQCD for D mesons, 
which have a different power counting than J/ip- 

For the simulations we use 278 quenched 12 3 x 
24 configurations at (3 = 5.7 generously provided 
by the UKQCD-Collaboration. To simulate the 
light quarks on the lattice we use the tadpole im- 
proved clover action with n c = 0.1434 from 
and k s \k = 0.1398 from ||. The Hamiltonian of 
the heavy quarks is expanded around its non rel- 
ativistic limit with tadpole improved gauge links 
and fields: 



H = H +5H, 
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This expansion is correct up to 0(Mq 2 ) with the 
last two terms being corrections for finite lattice 
spacing. In the quenched approximation heavy or 
light quark physics deliver different lattice spac- 
ings: from m p one gets a -1 = 1.103(11)(50) GeV 
§ and from the T-spectrum 1.4(1) GeV §J. 



2. MESON SPECTROSCOPY 

So far we have analysed correlators with the 
1 Sq and 3 Si quantum numbers and use two dif- 
ferent smearing functions for each value of Mq. 
When extracting masses we fit with two expo- 
nentials in case of the respective ground state and 
three in case of radial excitations. The position 
of the c and the b quark in the Mq range is fixed 
from the requirement that the D and B masses 
must coincide with their physical values in GeV 
H . To correct the meson mass for the omitted Mq 
in (|l|) we use perturbative shifts ||, which are in 
very good agreement with the non perturbative 
shifts from J/ij) and T 

In fig. |l| we give our results for the B spectrum 
together with preliminary results at (3 = 6.0 with 
the same action. For the comparison we set the 
scale from m p , since the uncertainty in how to de- 
termine a cancels out between both lattice results. 
The error bars represent the statistical uncertain- 
ties and those arising in the determination of the 
bare b quark mass. The horizontal lines give the 
experimental values, quoted from Q , with excep- 
tion of the radially excited B(2S), where we quote 
a preliminary DELPHI result || . Comparing the 
results, no significant scaling violation can be ob- 
served within the error bars of ks 10% in case of 
the hypcrfinc and B s — B splitting. 



2 



to 



1 — 
0.0 



I 600 
□ 

400 

m 

m 200 

5" 



o 1 — 

0.0 



100 



50 



1 — 
0.0 



m 20 



0.5 

a in GeV~ 



Figure 1. Scaling of the B spectrum with a. 



Our result for the D spectrum is given in fig. 
The arrows indicate the shift of the central val- 
ues when taking the average of the above quoted 
a- values instead of a(m p ), the other symbols are 
similar to fig. [I]. When comparing to the experi- 
mental values, the splittings between the strange 
and non strange mesons are in very good agree- 
ment and also the result for the radial excited 
B(2S) agrees within its fairly large error bar. 
However all hyperfine splittings are much too low 
when using a from m p . When using a larger 
value of a , for example from T, the hyper- 
fine increases but the agreement of the B s — B 
and D s — D is lost. The cause of this trouble 
might be a systematic error of the quenched ap- 
proximation, a quite large 1-loop correction to the 
coefficient in front of the a B term or neglected 
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Figure 2. Splittings in the D meson spectrum. 
The circles give the result for a' 1 = 1.10 GeV 
and the arrows for 1.25 GeV. The experimental 
values are given by horizontal lines. 



0(M~ 3 ) terms ineq. (0). This question definitely 
requires further investigation. 

It is interesting to have a look at the ratio of 
the hyperfine splittings. We find: 

D* - D , , D* — D s 

3-1(7), 15^-^ = 3.1(7), (4) 



B* -B 



B* - B s 



which has to be compared to an experimental 
value of 2.95(2) resp. 3.06(17). The error in (|) 
is dominated by the uncertainty of the perturba- 
tive mass shifts in the determination of the bare c 
quark mass. However this might be too conserva- 
tive and the agreement of the central values itself 
is quite impressive. 

3. DECAY CONSTANT 

In the effective theory one has to measure the 
matrix elements of the currents 
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in order to calculate the pseudoscalar decay con- 
stant up to 0(Mq 1 ). On the lattice for zero mo- 
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Figure 4. Scaling of /g s ■ The error bars contain 
statistical errors and the uncertainty of q* . 



Figure 3. Pseudoscalar decay constant. The 
connecting lines are to guide the eye. The po- 
sition of the B and the D meson are indicated. 
Statistical errors only are shown. 



mentum we have J 
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sation and matching to full QCD have recently 
been calculated to 1-loop order |||l(J. However 
as yet the matching scale q* has not been calcu- 
lated and we will give results for aq* = 1 and it. 
We will also apply the 0(aa s ) improvement term 
as introduced in the references. For the log(aAf) 
terms in the coefficients of and Jjj 1 we insert 
^log(aMo), which will give the physical decay 
constant for B and D mesons correctly. 

So far we have analysed the currents for k — 
0.1400 which is close to k s \k- From our expe- 
rience at smaller a we expect the effect of this 
mismatch to be 1% on /ps, which is negligible 
compared to the other errors. Our result is given 
in fig. ||. In the B region fpg is reduced quite 
largely by 15-25% due to renormalisation and im- 
provement. We obtain 



f Bs =211(7)(15)(37) MeV. 



(0) 



In parenthesis we give the statistical uncertainty, 
the uncertainty of q* and the shift of the central 
value for a" 1 = 1.25 GeV instead of 1.10 GeV. 
In fig. H we discuss the scaling of this quantity. 
Since the correct q* values need not to be equal 
at (3 = 5.7 and 6.0, we average over aq* = 1 and 
7r. The difference is included into the error bars. 



At tree level we observe a scaling violation of « 
30% which is reduced in case of the renormalised 
and improved result, so that they agree within 
error bars. Note that the result at larger a is 
more sensitive to q*. 

4. CONCLUSION 

We observe no significant scaling violations for 
the B spectrum and the renormalised fg s . Our 
results on the D spectrum look encouraging. 
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